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Abstract 

The results from the COBE*satenite are in close quantitative agreement with 
the predictions of the standard hot Big Bang model, suggesting that the 
Universe was once hot, dense and isothermal, giving a background radiation 
spectrum that is close to a perfect blackbody. The need to preserve this 
blackbody spectrum to match the spectrum observed by the FIRAS instrument 
on COBE places strong limits on events and scenarios occurring later than 
1 year after the Big Bang. The observation of intrinsic anisotropy of 
the microwave background by the DMR instrument on COBE provides a 
measurement of the magnitude of the gravitational potential fluctuations that 
existed 10^'^ years after the Big Bang. These potential perturbations were 
either produced in an inflationary epoch sometime during the first picosecond 
after the Big Bang, or else they are initial conditions dating from t = 0. The 
angular power spectrum of the correlations is compatible with the inflationary 
scenario. The observed amplitude of AT can be used to restrict theories of the 
fundamental structure of matter on micro-physical scales and at energies higher 
than that of the most energetic cosmic rays. The gravitational forces implied 
by the observed potential can produce the observed clustering of galaxies, but 
only if the density of the Universe is dominated by some form of transparent 
"dark" matter, which interacts weakly with photons. While the FIRAS and 
DMR instruments have 7° beams, the third instrument on COBE, DIRBE, has 
collected data with an 0.7° beam that will be used to improve our knowledge of 
the cosmic infrared background, but only after the strong foreground emissions 
from the interplanetary dust cloud and galactic dust are modeled and removed. 



1. Introduction 

The observables of modern cosmology include the Hubble expansion of the 
Universe; the ages of stars and clusters; the distribution and streaming motions 

* The National Aeronautics and Space Administration/Goddard Space Flight Center (NASA/GSFC) 
is responsible for the design, development, and operation of the Cosmic Background Explorer 
(COBE). Scientific guidance is provided by the COBE Science Working Group. GSFC is also 
responsible for the development of the analysis software and for the production of the mission data 
sets. 
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of galaxies; the content of the Universe (its mass density, composition, and the 
abundances of the hght elements); the existence, spectrum and anisotropy of the 
cosmic microwave background radiation; and other potential backgrounds in the 
infrared, ultraviolet, x-ray, and gamma-ray spectral regions. The purpose of the 
COBE mission is to make definitive measurements of two of these observable 
cosmological fossils: the cosmic microwave background (CMB) radiation and the 
cosmic infrared background (CIB) radiation. Since the discovery of the CMB^, 
many experiments have been performed to measure the CMB spectrum and spatial 
anisotropies over a wide range of wavelengths and angular scales. Fewer attempts 
have been made to conduct a sensitive search for a CIB radiation, expected to result 
from the cumulative emissions of luminous objects formed after the universe cooled 
sufficiently to permit the first stars and galaxies to form. 

The three scientific instruments on COBE are the Far Infrared Absolute 
Spectrophotometer (FIRAS), the Differential Microwave Radiometers (DMR), and 
the Diffuse Infrared Background Experiment (DIRBE). The purpose of FIRAS is to 
make a precision measurement of the spectrum of the CMB from 1 cm to 100 /xm. The 
purpose of DMR is to search for variations in the temperature of the CMB on angular 
scales larger than 7° at frequencies of 31.5, 53, and 90 GHz. The purpose of DIRBE is 
to search for a CIB by making absolute brightness measurements of the diffuse infrared 
radiation in 10 photometric bands from 1 to 240 fim and polarimetric measurements 
from 1 to 3.5 /i.m. The FIRAS and DIRBE instruments are located inside a 650 liter 
superfluid hquid helium dewar. See Boggess et al? for a full description of the COBE 
mission. 

The COBE mission is the product of many years of work by a large team 
of scientists and engineers. Credit for all of the results presented in this paper 
must be shared with the other members of the the COBE Science Working Group: 
Chuck Bennett, Ed Cheng, Eli Dwek, Mike Hauser, Tom Kelsall, John Mather, 
Harvey Moseley, Nancy Boggess, Rick Shafer, Bob Silverberg, George Smoot, Steve 
Meyer, Rai Weiss, Sam Gulkis, Mike Janssen, Dave Wilkinson, Phil Lubin, and Tom 
Murdock. Bennett et ai.^ is a review of the history of the COBE project and its 
results up to the discovery of the anisotropy by the DMR, but not including the 
latest FIRAS limits on distortions. 

Some of the members of this team have been working on COBE since the release 
of the NASA Announcement of Oppurtunities AO-6 and AO-7 in 1974, when the 
proposals for what became the COBE project were submitted. I have been working 
on COBE since the beginning of 1978. COBE was successfully launched at 14:34 UT 
on 18 November 1989 from Vandenberg AFB, California, and returned high quality 
scientific data from all three instruments for ten months until the liquid helium ran 
out. But about 50% of the instruments do not require liquid helium, and are still 
returning valuable scientific data in July 1993. While the COBE mission has been 
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successful, making the "discovery of the century" , one must remember that this work 
is based on the earher work (in the 20*^ century!) of Hubble'^ and Penzias and Wilson^, 
who discovered the expansion of the Universe and the microwave background itself. 
As a consequence of these two discoveries, one knows that the early Universe was 
hot and dense. When the density and temperature are high, the photon creation 
and destruction rates arc high, and are sufficient to guarantee the formation of a 
good blackbody spectrum. Later, as the Universe expands and cools, the photon 
creation and destruction rates become much slower than the expansion rate of the 
Universe, which allow distortions of the spectrum to survive. In the standard model 
the time from which distortions could survive is 1 year after the Big Bang at a redshift 
z ^ 10^-^. 

However, the expansion of the Universe will not produce a distortion by itself, since 
adiabatically expanding a blackbody radiation field results in another blackbody with 
a lower temperature. Thus the existence of a distorted spectrum in the hot Big Bang 
model requires the existence at time later than 1 year after the Big Bang of both an 
energy source and an emission mechanism that can produce photons that are now 
in the millimeter spectral range. Conversely, a lack of distortions can be used to 
place limits on any such energy source, such as decaying neutrinos, dissipation of 
turbulence, etc. 

At a time 10^-^ years^ after the Big Bang, a,t z ^ 1360, the temperature has 
fallen to the point where helium and then hydrogen have 50% (re) combined into 
transparent gases^. The surface of last scattering (9r/91n(l + z) = 1) occurs later, 
at z ^ 1160 or 10^'^ years after the Big Bang. The electron scattering which had 
impeded the free motion of the CMB photons until this epoch is removed, and the 
photons stream across the Universe. Before recombination, the radiation field at any 
point was constrained to be nearly isotropic because the rapid scattering scrambled 
the directions of photons. The radiation field was not required to be homogeneous, 
because the photons remained approximately fixed in comoving co-ordinates. After 
recombination, the free streaming of the photons has the effect of averaging the 
intensity of the microwave background over a region with a size equal to the horizon 
size. Thus after recombination any inhomogeneity in the microwave background 
spectrum is smoothed out. Note that this inhomogeneity is not lost: instead, it is 
converted into anisotropy. When we study the isotropy of the microwave background, 
we are looking back to the surface of last scattering ^ 10^'^ years after the Big Bang. 
But the hot spots and cold spots we are studying existed as inhomogeneities in the 
Universe before recombination. Since the 7° beam used by the DMR instrument 
on COBE is larger than the horizon size at recombination, these inhomogeneities 
cannot be constructed in a causal fashion during the epoch before recombination in 
the standard Big Bang model. Instead, they must be installed "just so" in the initial 

t With Ho = 50, n = 1. 
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conditions. In the inflationary scenario^ these large scale structures were once smaller 
than the horizon size during the inflationary epoch, but grew to be much larger than 
the horizon. Causal physics acting 10~^^ seconds after the Big Bang can produce the 
large-scale inhomogeneities studied by the DMR. 

2. COBE Orbit and Attitude 

COBE was launched at dawn on 18 November 1989 into a sun-synchronous 
orbit with an inclination of 99° and an altitude of 900 km. By choosing a suitable 
combination of inclination and altitude, the precession rate of the orbit can be set 
to follow the motion of the Sun around the sky at 1 cycle per year. Thus the line 
perpendicular to the orbit plane always points approximately toward the Sun. Since 
COBE was launched toward the South in the morning, the end of the orbit normal 
which points toward the Sun is at declination —9°. At the June solstice, when the 
declination of the Sun is +23°, there is a 32° deviation from the ideal situation where 
the Earth-Sun line is exactly perpendicular to the orbit. Since the depression of 
the horizon given by cose? = -R©/a, where a is the radius of the orbit of COBE, is 
29°, there are eclipses when COBE passes over the South pole for a period of two 
months centered on the June solstice. During the same season at the North pole the 
angle between the Sun and anti-Sunward hmb of the Earth is only 177°. It is thus 
not possible to keep both Sunlight and Earthshine out of the shaded cavity around 
the dewar. Since the Sun produces a much greater bolometric intensity than a thin 
crescent of Earth limb, the choice to keep the Sun below the plane of the shade is an 
obvious one. 

The pointing of COBE is normally set to maintain the angle between the spin 
axis and the Sun at 94°. This leaves the Sun 4° below the plane of the shade. During 
eclipse seasons this margin is shaved to 2° to reduce the amount of Earthshine into 
the shaded cavity where the instruments sit. The angle between the plane of the 
shade and the Sun is known as the "roll" angle. The "pitch" angle describes the 
rotation of the spin axis around the Earth-Sun line. This rotation is programmed for 
a constant angular rate of 1 cycle per orbit. The nadir is set to be close to the minus 
spin axis, but because the Sun can be up to 32° from the orbit normal and the Sun 
is kept 4° below the plane of the shade, the spin axis can be up to 36° away from the 
zenith. The uniform rate in pitch allows the nadir to oscillate by up to ±6° relative 
to the spin-Sun plane with a period of 2 cycles per orbit. This motion is similar to 
the effect of the inclination of the ecliptic on the equation of time that can be seen 
in an analemma. To keep the "wind" caused by the orbital motion of COBE from 
impinging on the cryogenic optics, the pitch is biased backward by 6°. The "yaw" 
motion of COBE is its spin. The spin period is about 73 seconds. 
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There are two types of instrument pointings on COBE. The DMR and DIRBE 
instruments have fields-of-view located 30° away from the spin axis. As a resuh they 
scan cycloidal paths through the sky, that cover the band from 64° to 124° from the 
Sun. This band contains 50% of the sky, and it is completely covered in a day by the 
DMR or DIRBE. This rapid coverage of the sky imphes a rapid scan rate: 2.5°/sec. 
The data rates are also high: the DIRBE transmits 8 rcadings/channel/sec, while 
the DMR transmits 2 readings/channel/sec. The FIRAS points straight out along 
the spin axis. As a result its field-of-view moves only at the orbital rate of 3.6°/min. 
Since FIRAS coadds interferograms over a collection period of about 30 to 60 seconds, 
this slow scanning of the sky is essential. 

Over the course of a year the Sun moves around the sky, and the scan paths of the 
instruments follow this motion. The leading edge of the wide scanned band for the 
DIRBE and DMR caught up to the trailing edge of the band at start of mission after 
4 months, leading to 100% sky coverage. For the FIRAS total sky coverage would in 
principle be obtained after 6 months, but because of gaps in the observations caused 
by calibration runs and lunar interference, the actual sky coverage for FIRAS, defined 
as having the center of the beam within a pixel, is about 90%. The coverage gaps are 
sufficiently narrow, however, that every direction on the sky was covered by at least 
the half power contour of the beam. 

The attitude determination system for COBE uses infrared star sightings from 
DIRBE to correct a gyro propagated attitude solution. This method requires 3 
working gyros that are not co-planar. On 10 July 1993, COBE suffered its fourth 
gyro failure out of six gyros, which will degrade the accuracy of the attitude solutions 
by about a factor of 5, to a la accuracy of about 0.1°. (One of the failed gyros is 
noisy but operable.) This accuracy will be sufficient for the DMR experiment with its 
7° beam, provided that the errors are not systematic. Current planning is to operate 
COBE until December 1993, which will complete four years of observations with the 
DMR and the short wavelength channels of the DIRBE. 

3. DIRBE Meeisurements 

The Diffuse Infrared Background experiment (DIRBE) on COBE is designed to 
search for a cosmic infrared background (GIB). Models of the GIB radiation take 
light from distant, unresolved galaxies, and redshift it into the near infrared bands'^; 
or absorb the starhght on dust, reradiate it in the far-infrared, and then redshift it 
into the sub- millimeter band^. These models predict fluxes that are approximately 
two orders of magnitude lower than the emission from the Milky Way, which makes 
detection of the GIB a difficult task, unless observations are made at wavelengths 
where the redshift carries the energy from distant galaxies into a region with little 



6 



Milky Way emission. The interplanetary dust cloud, which scatters sunlight to give 
the zodiacal light, absorbs much more power than it scatters and reradiates this 
power in the mid-infrared. Thus from 5 /xm to 100 /xm the CIB is overwhelmed when 
observing from a distance 1 au from the Sun. However, there are windows between 
the scattered and thermal zodiacal light at 3.5 fim and between the interstellar dust 
in the Milky Way and the CMBR at 200-300 /im where the CIB may be detectable. 

The DIRBE instrument is an off-axis Gregorian telescope with a 19 cm diameter 
primary. The Gregorian design has an intermediate image of the sky, which allows 
for the placement of a stray-light stop between the primary and the secondary, and 
the secondary forms an image of the primary where a second stray-light stop is 
located. Baffle vanes in the telescope tube absorb radiation coming from far off-axis, 
so the main path for stray light involves scattering on the primary mirror when it is 
illuminated by a slight off-axis source. Great care was taken to avoid contamination 
and to preserve the cleanliness of the optics, so the off-axis response of the DIRBE is 
very small. The DIRBE beam is an 0.7° x 0.7° square. Dichroics and aperture division 
are used to gives simultaneous measurements in 10 spectral bands in the same field of 
view. The 10 bands in DIRBE are the ground-based photometry bands J[1.2], K[2.3], 
L[3.4], and M[4.9]; the four IRAS bands at 12, 25, 60, and 100 /im; and two long 
wavelength bands at 140 and 240 /xm. The shortest wavelengths (JKL) bands, which 
observe scattered light from the interplanetary dust cloud, have 3 detectors per band 
with polarizers to measure linear polarization. A tuning fork running at 32 Hz serves 
as a chopper for the DIRBE. This gives a well-determined zero level by comparing 
the sky signal to the 2 K interior of the DIRBE, allowing the search for an isotropic 
CIB to proceed. A shutter blocks the optical path less frequently, so the accuracy 
of the zero level can be checked. When the shutter is in the beam, its back side is 
a mirror which refiect light from a commandable internal reference source (IRS) into 
the detectors, providing a gain calibration several times per orbit. During the 10 
months of cryogenic operation, the achieved instrument rms sensitivity per field of 
view was A/(A) = (1.0, 0.9, 0.6, 0.5, 0.3, 0.4, 0.4, 0.1, 11.0, 4.0) x 10"^ W m'^ 
sr~^, respectively for the ten wavelength bands listed above. These noise levels are 
considerably smaller than the total fiuxes measured by DIRBE at the South Ecliptic 
Pole given in Table |I[ Note that the predicted extragalactic fiux is 3 x 10~^ Wm~^sr~^ 
at both 3.5 /xm^ and 240 /im^ which is 20% of the total SEP fiux at 3.5 /xm and 43% 
at 240 yum. 

The final determination of the magnitude of the extragalactic background will 
depend on modeling and removing the zodiacal light. The DIRBE beam is aligned 
30° from the spin axis of COBE, so it swings from 64° elongation to 124° elongation in 
every rotation of the satellite. In principle, the variation of intensity with elongation 
can be used to determine the magnitude of the zodiacal emission, just as the variation 
of atmospheric extinction with secant of the zenith angle allows one to determine the 
fiux hitting the top of the atmosphere. However, the shape of the atmosphere is well 
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Reference 


A 


A/a 


^^(A) 




(um) 


(10-^ W m-2 sr-i) 


K 


DIRBE 


1.2 


8.3 ± 3.3 


351 


{South Ecliptic Pole) 


2.3 


3.5 ± 1.4 


216 




3.4 


1.5 ± 0.6 


140 




4.9 


3.7 ± 1.5 


109 




12. 


29. =fc 12 


56 




22. 


21. ± 8 


34 




55. 


2.3 ± 1 


15 




96. 


1.2 ± 0.5 


9.2 




151. 


1.3 ± 0.7 


6.6 




241. 


0.7 ± 0.4 


4.5 



Table 1: Brightness and Planck Brightness Temperature of the Diffuse Infrared Sky 

known, while the shape of the zodiacal dust cloud is poorly determined. If the cloud 
is in equiUbrium, then Jeans' theorem states that the phase space density should be a 
function of the integrals of the motion: total angular momentum, L, the z component 
of the angular momentum, L^, and the energy, E. The total angular momentum and 
energy define the radius and eccentricity of an orbit. If the eccentricity distribution 
can be ignored, one is left with a dependence on radius. The ratio of to L defines 
the inclination of the orbit. Thus the cloud should be describable by a function 
f{r,i), which suggests that fan- shaped models should be used. While f{r,i) allows 
for the "bands" seen by IRAS, the comet trails are clearly non-equilibrium structures 
outside the scope of these models. Moreover, the center from which r is measured 
and the axis from which i is measured depend on planetary perturbations^ and are 
both functions of r. Clearly a lot of effort will be needed to achieve the hoped for 
removal of the zodiacal light to an accuracy of 1% of the peak emission. 

After the depletion of the liquid helium, the DIRBE JKLM channels continued to 

operate. The interior temperature of the Dewar appears to have stabilized at ~ 55 K, 
which is adequate for the InSb detectors used in these bands. While the "blanked-off" 
noise in these channels has risen due to the increased temperature, the noise while 
observing the sky is dominated by the confusion noise due to the multitude of faint 
stars. These "warm" data have been used to provide star sightings for the attitude 
determination. 

4. DMR Observations 

The Differential Microwave Radiometers (DMR) experiment on COBE is designed 
to measure small temperature differences from place to place on the sky. The DMR 
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consists of three separate units, one for each of the three frequencies of 31.5, 53 and 90 
GHz. The field of view of each unit consists of two beams that are separated by a 60° 
angle that is bisected by the spin axis. Each beam has a 7° FWHM. The DMR is only 
sensitive to the brightness difference between these two beams. This differencing is 
performed by a ferrite waveguide switch that connects the receiver input to one horn 
and then the other at a rate of 100 Hz. The signal then goes through a mixer, an IF 
amplifier and a video detector. The output of the video detector is demodulated by a 
lock-in amplifier synchronized to the input switch. The difference signal that results 
is telemetered to the ground every 0.5 seconds. Each radiometer has two channels: A 
and B. In the case of the 31.5 GHz radiometer, the two channels use a single pair of 
horns in opposite senses of circular polarization. In the 53 and 90 GHz radiometers 
there are 4 horns, and all observe the same sense of linear polarization. 

The basic problem in the DMR data analysis is to construct a map of the sky using 
the 6 X 10'' differences that are collected each year in each of the six channels. There 
are only about 10^ beam areas on the sky, so the problem is highly over-determined. 
We have chosen to analyze the data using 6,144 pixels to cover the sky. These are 
approximately equal area, and arranged in a square grid of 32 x 32 pixels on each of 
the 6 faces of a cube. Within each pixel we assume that the temperature is constant, 
so we are modeling the sky with a staircase function. Therefore the basic problem can 
be represented as a least-squares problem with 6 x 10 equations in 6,144 variables. 
One way to solve this problem is to average all of the observations of each pixel. This 
gives the temperature of the pixel relative to its "reference ring" of pixels 60° away. 
Figure |I| shows the reference ring for pixel 4342 which contains the galactic center. 
Then correct the average for each pixel by the weighted mean of the temperatures in 
the reference ring. The correction must be iterated using the corrected temperatures 
until the process converges, but the convergence is slow. In practice, the DMR maps 
are computed by inverting the 6144 x 6144 matrix of normal equations. Fortunately 
this matrix Aij is sparse and symmetric. Thus only 9 x 10^ elements need to be kept. 
Figure |I| can be viewed as a map made by plotting Aij at the position of the i*^ pixel 
while j is fixed at j = 4342. The diagonal An shown in Figure |^ shows the number 
of observations per pixel. This figure is plotted in galactic coordinates to better show 
the structure of the heavily covered rings 26° and 34° from the ecliptic poles. The 
coverage near the North ecliptic pole is reduced because the Earth limb rises above 
the sun shade over the North pole during the May- July eclipse season. 



5. DMR Results 



The basic DMR result is the discovery of an intrinsic anisotropy^'' of the microwave 
background, beyond the dipole anisotropy^^'^^. This anisotropy, when a monopole and 
dipole fit to \b\ > 20° is removed from the map, and the map is then smoothed 
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Row 4342 of Matrix A 



-125.0 -100.0 -75.0 -50.0 -25.0 0.0 

Figure 1: The column through pixel 4342 of the matrix A for the 53 GHz A channel. 
The central pixel contains the galactic center, and is saturated with a value of 10,599. 
Ecliptic coordinates are used for this plot. 



to a resolution of ~ 10°, is 30 /iK. The correlation function of this anisotropy 
is well fit by the expected correlation function for the Harrison-Zeldovich^^'^^'^^ 
spectrum of primordial density perturbations predicted by the inflationary scenario, 
and the amplitude is consistent with many models of structure formation^^. Galactic 
emission^'', systematic errors^^ and known extra-galactic sources^^ cannot explain the 
observed signal. 

Even though the noise per pixel is much larger than the cosmic anisotropy, the 
detection of anisotropy has a high statistical confidence level because of the large 
number of pixels in the DMR maps. To illustrate this point. Figure |^ shows the 
unsmoothed A and B side 53 GHz maps from the initial data products, plus the 
smoothed difference and sum maps. Compare these to a test case with a known 
map, degraded to a comparable signal to noise ratio, shown in Figure |. Clearly 
a pattern of some sort is starting to emerge in the smoothed sum map. Figure 
I shows the two auto-correlations and the cross-correlation that can be computed 
from the unsmoothed test case maps. The solid curve in this figure is the known 
correlation function based on the known input map. (A noiseless high resolution 



picture of the input map is in Figure Il3l) Clearly the DMR maps show statistically 
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Coverage Map for 53 A Channel 




0. 10000. 20000. 30000. 40000. 50000. 



Figure 2: The diagonal of the matrix A for the 53 A channel, otherwise known as the 
coverage map, in galactic coordinates. 

significant correlated structure. The dashed curve is the correlation function of a 
Harrison-Zeldovich model with equal power on all scales, smoothed to the DMR 
resolution, with the amplitude adjusted to match the test case. The relatively good 
agreement between the dashed curve and the test case correlation function shows 
that the good agreement between the DMR correlation function and the Harrison- 
Zeldovich model^°, while supporting inflation, does not prove inflation. 

The inflationary scenario also predicts that the temperature fluctuations should 
have a Gaussian distribution. This can be tested by looking at higher-order moments 
of the maps such as three point correlation functions, or by directly studying the 
probability distributions. The random measurement noise from the radiometers is still 
sufficient to interfere with these studies to a great extent, but a preliminary analysis 
finds that the Gaussian model is consistent with the observations. Figure |^ shows 
histograms constructed from the smoothed DMR maps. To minimize the radiometer 
noise, the thermodynamic AT's derived from the 53 and 90 GHz channels have been 
averaged together with weights of 0.6 and 0.4 respectively. The histogram of the 
{A—B)/2 difference map is not expected to be Gaussian, because of the large variation 
in integration time per pixel seen in Figure |^. However, a Gaussian convolved with 
the {A — B)/2 histogram should match the histogram of the {A + B)/2 map if the sky 
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53 GHz A 53 GHz B 




-1.250 -0.250 0.750 1.750 2.750 3.750 -1.250 -0.250 0.750 1.750 2.750 3.750 



Smoothed (A-B)/2 



Smoothed (A+B)/2 




-.3000 -.1000 0.1000 0.3000 0.5000 0.7000 -.3000 -.1000 0.1000 0.3000 0.5000 0.7000 



Figure 3: The unsmoothed 53 GHz A & B maps from the initial COBE data products, 
and the smoothed sum and difference maps. 



has a Gaussian distribution. As Figure § shows, the "Gaussian" curve, which is this 
convolution, matches the observed histogram well. Detailed calculations show that 
models which predict non-Gaussian temperature fluctuations, such as cosmic strings^^ 
or global monopoles^^, can also match the observed histogram. This convergence 
toward Gaussian probability distributions is caused by the 10° effective resolution 
of the smoothed maps, which averages many small patches together. Thus cosmic 
strings, which predict dramatically non-Gaussian behavior (sharp edges) on small 
scales, give almost Gaussian results after smoothing. However, some models could be 
rejected: Figure ^ shows the histograms from the test case maps in Figure ^. Here 
the non-Gaussian structure is correlated over large distances, so the smoothing does 
not erase it. 

The ratio of the gravitational potential fluctuations seen by the DMR through 
the Sachs-Wolfe^^ effect to the gravitational potentials inferred from the bulk flows of 
galaxies using the POTENT method^^, shows that the Harrison-Zeldovich spectrum 
has the correct slope to connect the DMR data at 3x 10^ km/sec scales to the bulk flow 
data at 6000 km/sec scale. This comparison of the horizon-scale perturbations seen 
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Test Case A Test Case B 




-1.250 -0.250 0.750 1.750 2.750 3.750 -1.250 -0.250 0.750 1.750 2.750 3.750 




Figure 4: The unsmoothed test case A & B maps made by degrading a known pattern 
to the DMR resolution and signal to noise ratio, and the smoothed sum and difference 
maps. 



by COBE and the "large-scale" {60/h Mpc, where h = iJo/(100 km s'^ Mpc"^)) 
structure data gives best determination of the slope of the power spectrum of the 
primordial perturbations, even though it is model-dependent. COBE data alone 
cover such a small range of scales that slope determinations are uncertain. 

Wright et al}^ discussed the implications of the COBE DMR data for models 
of structure formation, and selected 4 models from a large collection^* for detailed 
discussion: a "CDM" model, with Ho — 50 km s~^ Mpc~^, 0,cdm — 0.9, and 

= 0.1; a mixed "CDM+HDM" model, with Ho = 50 km s"^ Mpc-\ (^cdm = 0.6, 
^HDM = 0.3 (a 7 eV neutrino), and Qb = 0.1; an open model, with Ho = 
100 km s-^ Mpc-\ flcDM = 0.18, and Qb = 0.02; and a vacuum dominated model, 
with Ho = 100 kms-^ Mpc"\ Qcdm = 0.18, Qb = 0.02, and Q„„c = 0.8. The 
vacuum dominated model and especially the open model have potential perturbations 
now that are too small to explain the POTENT bulk flows^^. Note that all of these 
models are normalized to the COBE anisotropy at large scales at z = 10'^, but in the 
open and vacuum-dominated models the potential perturbations go down even at the 
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Figure 5: The auto-correlations of the test case A & B maps (sohd dots), the cross- 
correlation (open circles), and the correlation function of the known input pattern 
at DMR resolution. Also shown as a dashed curve is the correlation function of an 
n = 1 Harrison- Zeldovich model. 



largest scales which are not affected by non-linearities or the transfer function. The 
small-scale anisotropics predicted by these models are close to the upper limits, and in 
the case of the 1.5° beam South Pole experiments^ (SP) they exceed the current data. 
Since the DMR datum is at £e// = 4 while the SP datum is at ie/f = 44, a tilted model 
which replaces the Harrison-Zeldovich n = 1 in P{k) oc A;" with an n ~ 0.5 would 
eliminate this discrepancy. But such a low n would destroy the agreement between 
the bulk flow velocities and the COBE {Q'rmsT'^- In fact the POTENT bulk flows 
are inconsistent with the SP AT data^^, while other bulk flow determinations are 
even larger^^. 

Wright^^ has computed the power spectrum of the DMR maps using the Hauser- 
Peebles method to allow for incomplete sky coverage caused by masking out the 
galactic plane region. Figure ^ shows the power in each multipole normalized to a 
pure n = 1 (Harrison-Zeldovich) power law with (Q'rhjs)^'^ ~ ^"^ f^^- Also shown 




Figure 6: The histogram of the smoothed 1 year 53+90 GHz (A-B)/2 map, its 
convolution by a Gaussian, which matches the observed histogram of the (A+B)/2 
map. 

are results from various small and medium scale AT experiments. For perturbations 
at scales smaller than the horizon at the end of radiation dominance, dynamical 
effects boost the expected anisotropy. At very small scales, the finite thickness of the 
recombination surface filters out most of the anisotropy. The models shown in Figure 
^ are an n = 0.85 tilted CDM model^^ (lower curve, with Ho = 50 km s~^ Mpc~^ 
and Qsh^ = 0.0125), and the same model with the n = 0.96 expected from chaotic 
0'^ inflation. The POTENT bulk flow, if translated into a AT/ in Figure ||, would 
correspond to a point on the higher model curve at £ ~ 100. The Lauer & Postman^^ 
bulk flow would correspond to a AT/ that is ~ 40 times higher than the higher model 
at £ ~ 40. Clearly the tilt of the primordial power spectrum cannot be deterined 
without much better measurements of small scale AT and bulk flow velocities. Note 
that the RING experiments^ is known to be contaminated by discrete sources, but the 
power spectrum of uncorrelated point sources is cx i"^ in Figure so this population 
has a negligible effect on the data at lower ts for sources which emit a spectrum 
similar to the spectrum of the Milky Way. When comparing these data to the COBE 




Figure 7: The histogram of the smoothed (A-B)/2 test case map, and its convolution 
by a Gaussian, which fails to match the observed histogram of the (A+B)/2 map. 

DMR data the corrections for the galactic plane cut and the non-Gaussian shape of 
the DMR beam should be used^°. 



6. FIRAS Meaisurements 



The Far InfraRed Absolute Spectrophotometer instrument on COBE is a 
polarizing Michclson^^ interferometer. The optical layout is symmetrical, and it has 
two inputs and two outputs. If the two inputs arc denoted SKY and ICAL, then the 
two outputs, which are denoted LEFT and RIGHT, are given symbolically as 

LEFT = SKY -ICAL (1) 
RIGHT = ICAL -SKY (2) 

The FIRAS has achieved its incredible sensitivity to small deviations from a blackbody 
spectrum by connecting the ICAL input to an internal calibrator, a reference 
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Figure 8: The angular power spectrum of the CMBR, normahzed to an {Q\msY'^ ~ 
17 fiK Harrison-Zeldovich model. Data from i = 2.. 13 are from COBE. Other data 
and upper limits are from left to right: Ulisse^^, Tenerife^"^, SP IP^, PYTHON^'^, 
MSAM35, two points from MAX^^'^^, MSAM again, RING=^« and OVRO^^ 
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blackbody that can be set to a temperature close to the temperature To of the sky. 
Thus this "absolute" spectrophotometer is so successful because it is differential. In 
addition, each output is further divided by a dichroic beamsplitter into a low frequency 
channel (2-21 cm"^) and a high frequency channel (23-95 cm~^). Thus there are four 
overall outputs. These are labeled LL (left low) through RH (right high). 

Each output has a large composite bolometer that senses the output power 
by absorbing the radiation, converting this power into heat, and then detecting 
the temperature rise of a substrate using a small and sensitive silicon resistance 
thermometer. In order to minimize the specific heat of the bolometer, which 
maximizes the temperature rise for a given input power, the substrate is made 
of the material with the highest known Debye temperature: diamond. Diamond 
is transparent to millimeter waves, however, so an absorbing layer is needed. 
Traditionally a layer of bismuth has been used as the absorber on composite 
bolometers, but the FIRAS detectors use a thin layer of gold alloyed with chromium. 
The absorbing layer needs to have a surface impedance about one-half of the 377 
Ohms/n impedance of free space in order to absorb efficiently, and this surface 
impedance can be obtained with a much thinner layer of gold rather than bismuth. 
Again, this serves to minimize the specific heat of the detector. In order to detect the 
longest waves observed by FIRAS, the diameter of the octagonal diamond substrate 
is large: about 8 mm. Each detector sits behind a compound parabolic concentrator, 
or Winston cone, which provides a full tt steradians in the incoming beam. As a 
result, the FIRAS instrument has a large etendue of 1.5 cm^ sr. 

Since the FIRAS is a Michelson interferometer, the spectral data are obtained in 
the form of interferograms. Thus the LEFT output is approximately 

Il{x) = cos(27ri/x)G'(i/) (^I, + J2 ei{v)B,{Ti) + u}j dv (3) 

The index i above runs over all the components in the FIRAS that had thermometers 
to measure Xi. These include the ICAL, the reference horn that connects to the 
ICAL, the sky horn, the bolometer housing, the optical structure of the FIRAS, and 
the dihedral mirrors that move to provide the variation in path length difference 
X. The ej's are the effective cmissivitics of the various components. The ICAL 
itself has e ?=a —1, while the sky and reference horns have e's of ± a few percent 
in the low frequency channels. The offset term was observed during flight to be 
approximately 10~^ exp(— t/r)i?,y(T[/) with a time constant r of two months and a 
temperature Tjy 15 K. The SKY input Ii, above can be either the sky or an external 
calibrator. The XCAL is a movable re-entrant absorber that can be inserted at the 
top of the sky horn. The combination of sky horn plus XCAL forms a cavity with an 
absorbtivity known to be > 0.9999 from measurements, and believed to be ^0.99999 
from calculations. With the XCAL inserted during periodic calibration runs, the 
SKY input is known to be B^{Tx)- By varying Tx and the other Tj's, the calibration 
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Figure 9: Co-added FIRAS interferograms from the Left Low channel taken during 
flight with the ICAL at 2.759 K while looking at the sky or the XCAL at 2.701, 2.733 
or 2.750 K. 



coefficients have been determined^". Figure ^ shows the interferograms from flight 
used to determine the preliminary FIRAS spectrum^^. 

Once the calibration coefficients are known, the sky data can be analyzed to 
determine b) the intensity of the sky as a function of frequency, galactic longitude 
and galactic latitude. This is a "data cube". The data from each direction on the 
sky can be written as a combination of cosmic plus galactic signals: 



T^(Z,6,oo) 



{B, iTo + AT{l,b)) + Ah] 



+ 



(4) 



where Tu{1, b, s) is the optical depth between the Solar system and the point at distance 
s in the direction (/, b) at frequency z/, AI^, is an isotropic cosmic distortion, and 
AT{l,b) is the variation of the background temperature around its mean value Tq. 
This equation can be simplified because the optical depth of the galactic dust emission 
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is always small in the millimeter and sub- millimeter bands covered by FIRAS. 

b) ^ B, (To + Ar(/, b)) + M, + j b, s)ds (5) 

Even in the optically thin limit, some restrictive assumptions about the galactic 
emissivity j^, are needed, since the galactic intensity J j,y{l,b, s)ds is a function of 
three variables, just like the observed data. The simplest reasonable model^^ for the 
galactic emission is 

''Ul,b,s)ds = G{l,b)g{u). (6) 



This model assumes that the shape of the galactic spectrum is independent of direction 
on the sky. It is reasonably successful except that the galactic center region is clearly 
hotter than the rest of the galaxy. The application of this model proceeds in two steps. 
The first step assumes that the cosmic distortions vanish, and that an approximation 
(7o(^) to the galactic spectrum is known. A least squares fit over the spectrum in each 
pixel then gives the maps AT(/, b) and G{1, b). The high frequency channel of FIRAS 
is used to derive G{1, b) because the galactic emission is strongest there. An alternative 
way to derive G{l,b) is to smooth the DIRBE map at 240 /im to the FIRAS 7° beam. 
This DIRBE method has been used in the latest FIRAS spectral results"'^'"'^'''^. The 
second step in the galactic fitting then derives spectra associated with the main 
components of the millimeter wave sky: the isotropic cosmic background, the dipole 
anisotropy, and the galactic emission. This fit is done by fitting all the pixels (except 
for the galactic center region with |6| < 20 and |Z| < 40) at each frequency to the 
form 

I,{l,b) = hiu) + D{iy) cose + G{l,b)g{i^). (7) 

The spectra of the anisotropic components derived from this fit applied to the entire 
mission data set are shown in Figures and 11. Figure 11 shows the residual^^ after 



the predicted dipole spectrum ^TdBy/dT evaluated at To is subtracted from -D(z/). 

But the FIRAS calibration model is so complicated and the cosmic distortions, if 
any, are so small, that there are still systematic uncertainties that limit the accuracy 
of the isotropic spectrum /o(z/) derived from the whole FIRAS data set. Therefore the 
best estimate of /o(i^) comes from the last six weeks of the mission, when a modified 
observing sequence consisting of 3.5 days of observing the sky with the ICAL set 
to null out the cosmic spectrum, followed by 3.5 days of observing the XCAL with 
its temperature set to match the sky temperature Tq. Six cycles of this alternation 
between sky and XCAL were obtained before the helium ran out. By processing 
both the sky data and the XCAL data through the same calibration model, and 
then subtracting the two spectra, almost all of the systematic errors cancel out. The 
residual in the high galactic latitude sky is then modeled as 

hisky, \b\ > 6,) - I^iXCAL) = AJ, + 6T—^ + Gg{u). (8) 

oT 
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The result of a least squares fit to minimize by adjusting 8T and G is shown 
in Figure [1^. Since all of the data taken during the last six weeks of the mission 
were taken in a single scan mode, the mechanical resonance frequency of the mirror 
carriage occurs at a fixed spectral frequency, leading to the large error bar on the 
point at 11 cm^^. The data points show the residuals while the curves show the 
"uninteresting" models dBy/dT and (7(2^), and the "interesting" y and // models. 
It is important to remember that any cosmic distortion with a spectral shape that 
matches one of the uninteresting curves will be hidden in this fit. The maximum 
residual between 2 and 20 cm~^ is 1 part in 3000 of the peak of the blackbody, and 
the weighted RMS residual is 1 part in 10^ of the peak of the blackbody. While 
these residuals are small, they are nonetheless more than twice the residuals that 



one would expect from detector noise alone. The error bars in Figure [ly have been 
inflated by a constant factor to give = 32 for the 32 degrees of freedom in the fit. 
Limits on cosmo logical distortions can now be set by adding terms to the above fit, 
finding the best fit value with minimum x^, and then finding the endpoints of the 
95% confidence interval where = xLin + 4- As an example of such a fit, suppose 
that the sky in fact was a gray-body with an emissivity e = 1 + e, where e is a small 
parameter. Since the XCAL is a good blackbody, this model predicts a residual of the 
form eByiTo). Because B^iTo) peaks at a lower frequency than dB^/dT, this model 
is sufficiently different from the "uninteresting" models and a tight limit on e can 
be set: |e| < 0.00041 (95% confidence). The result of fits for the usual cosmological 
suspects are \y\ < 2.5 x 10~^ and < 3.3 x 10""^ with 95% confidence. 

The absolute temperature of the cosmic background. To, can be determined two 
ways using FIRAS. The first way is to use the readings of the germanium resistance 
thermometers in the XCAL when the XCAL temperature is set to match the sky. 
This gives To = 2.730 K. The second way is to measure the frequency of the peak 
of dBy/dT by varying Tx a small amount around the temperature which matches 
the sky, and then apply the Wien displacement law to convert this frequency into a 
temperature. This calculation is done automatically by the calibration software, and 
it gives a value of 2.722 K for To. Three additional determinations of To depend on the 
dipole anisotropy. For either FIRAS or DMR, the spectrum of the dipole anisotropy 
can be fit to the form 

c oT 

Since the velocity of the solar system with respect to the CMB is not known a priori, 
only the shape and not the amplitude of the dipole spectrum can be used to determine 
To. For FIRAS, this analysis gives To = 2.714 ± 0.022 while for DMR it gives 
To = 2.76 ±0.18 K^^. The DMR data analysis keeps track of the changes in the dipole 
caused by the variation of the Earth's velocity around the Sun during the year. In 
this case the velocity v is known, so To can be determined from the amplitude of 
the change in the dipole, giving To = 2.75 ± 0.05 K^^. The final adopted value is 
2.726 ± 0.010 K"^^ (95% confidence), which just splits the difference between the two 
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Figure 10: FIRAS residuals (points) after fitting for a blackbody and galactic 
emission, with an 0.5 mK AT, a A(csc|fe|) = 0.25 galaxy curve, y — 2.5 x 10~^, 
and fjL — 3.3 X 10~^ shown for comparison. Full-scale is 0.1% of the peak of the CMB. 

methods based on the FIRAS spectra. The dipole-based determinations of To are less 
precise but provide a useful confirmation of the spectral data. 



7. FIRAS Interpretation 



When interpreting the implications of the FIRAS spectrum on effects that could 
distort the microwave background, it is important to remember some basic order-of- 
magnitude facts. Theoretical analyses of light element abundances gives limits on the 
current density of baryons through models of Big Bang NucleoSynthesis (BBNS). The 
ratio of baryons to photons^^ is r/ = (3.5 ±0.7) x 10"^*^, while another determination^^ 
gives f] = (3.95 ± 0.25) x 10"^'^. As a result, the effect of Lyman a photons from 
recombination is negligible, simply because there are so few hydrogen atoms per 
photon. Since the number density of photons is well determined by the FIRAS 
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Figure 11: FIRAS Dipole spectrum residuals from an amplitude of 3.343 mK applied 
to To = 2.726 K (points). Full-scale is 0.04% of the peak of the CMB. The curve 
shows the best fit dipole color temperature of 2.714 K with an amplitude of 3.379 
mK. 

spectrum 

AT = 87rC(3)r(3) j =410±5cm-3 (10) 

the baryon density is well determined: 1.62 x 10~^ cm~^. This gives n^/i^ = 0.0144. 
Not all of these baryons are protons, however. For a primordial helium abundance 
of 23.5% by mass, one has 13 protons and 1 a particle in 17 baryons. This gives an 
electron density of rie^ = 1.43 x 10^^ cm~^. For a fully ionized Universe, the optical 
depth for electron scattering over a path length equal to the Hubble radius c/Ho is 
0.00088//1. 

The expansion of the Universe H varies with time and redshift in a way that 
depends on ratios of the current matter, radiation and vacuum densities to the critical 
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density pc — SH'^/SnG. Let these ratios be Q^, Q,., and Q^. 



Then 



H{z) = (1 + z)Ho^{l -ilm-^r- ^v) + ^m{l + z) + 0,(1 + zf + ^^{1 + Z^ 

(11) 

This equation is not exact for massive neutrinos which shift from being part of 
into part of fi,. as z increases. The value of Q^/i^ = (4.16 ± 0.06) x 10~^ is fairly 

well-determined for the case of three massless neutrinos. Thus if 0^ ~ 1 now, 
the Universe becomes radiation dominated for z > z, 



eq 



2Ah^ X 10"^ when the 
background temperature was kTf.q = 5.65/i^ eV. In the radiation dominated epoch, 
H = 0.645(1 + zf km/sec/Mpc or // = 2.09(1 + zf x lO-^o g-^ 

For early epochs when the Universe was ionized, electron scattering is the 
dominant mechanism for transferring energy between the radiation field and the 
matter. As long as the electron temperature is less than about 10^ K, the effect 
of electron scattering on the spectrum can be calculated using the Kompaneets 
equation' 



48. 



dn 
dy 



— X 



,d_ 
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4 / 2 

X \n + n + -g^ j 



where n is the number of photons per mode (n 
X — hu/kTe, and the Kompaneets y is defined by 



dy 



kTe 



UeUTCdt. 



(12) 

l/(e^ — 1) for a blackbody), 

(13) 



Thus y is the electron scattering optical depth times the electron temperature in units 
of the electron rest mass. Note that the electrons are assumed to follow a Maxwellian 
distribution, but that the photon spectrum is completely arbitrary. 

Since the Kompaneets equation is describing electron scattering, which preserves 
the number of photons, one finds that the y derivative of the photon density N 
vanishes: 



dN 
dy 



oc 



x^—dx 
dy 



-I 



= 



d_ 
dx 



A ( 2 "^nX 
X I n + n + ^ I 



dx 



(14) 



The stationary solutions of the equation dn/dy = are the photon distributions 
in thermal equilibrium with the electrons. Since photons are conserved, the photon 
number density does not have to agree with the photon number density in a blackbody 
at the electron temperature. Thus a more general Bose-Einstein thermal distribution 
is allowed: n — l/(exp(x + //) — 1). This gives dn/dy — for all Since the 
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Bose-Einstein spectrum is a stationary point of the Kompaneets equation, it is the 
expected form for distortions produced at epochs when 

(1 + = aTne,o^^(l + zf>l (15) 

For the fi^/i^ given by BBNS, this redshift Zy where this inequahty is crossed is well 
within the radiation dominated era, with a value Zy = lO^ Vv^TOTV?. 

There is a simple solution to the Kompaneets equation with non-zero dn/dy which 
gives the Sunyaev-Zeldovich^^ or y spectral distortion. This simple case occurs when 
the initial photon field is a blackbody with a temperature which is below the 
electron temperature. Letting / = Tg/T^, we find that the initial photon field is 
given by n = l/(exp(/x) — 1). Therefore 

Thus the distortion has a Sunyaev-Zeldovich shape but is reduced in magnitude by a 
factor (1 — Ty/Tg). Defining the "distorting" y as 

dyo = 5 — neaxcdt (17) 

we find that the final spectrum is given by a frequency-dependent temperature given 
by 

T„ = n l + yol -41+.. 



where x = hv/kTo- The FIRAS spectrum in Figure |T0| shows that < 2.5 x 10 



The energy density transferred from the hotter electrons to the cooler photons in 
the y distortion is easily computed. The energy density is given by f/ oc / x^ndx so 

'xl{x^^^\dx (19) 



dyo J dx \ dx 
which when integrated by parts twice gives 

dU f f ^dn\ 



dyo J \ ox J 

x^ndx = 4U. (20) 



Thus the limit on yD gives a corresponding limit on energy transfer: AU/U < 10 ^. 
Any energy which is transferred into the electrons at redshift s z > 1 where the 
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Compton cooling time is less than the Hubble time will be transferred into the 
photon field and produce a y distortion. Since there are 10^ times more photons 
than any other particles except for the neutrinos, the specific heat of the photon gas 
is overwhelmingly dominant, and the electrons rapidly cool (in a Compton cooling 
time) back into equilibrium with the photons. The energy gained by the photons 
is Af/ = 4yaT^ which must be equal to the energy lost by the electrons and ions: 
1.5(ne + ni)kATe. Since y = axneikTe/ mec'^)cAt we find the Compton cooling time 



for an ion to electron ratio of 14/15. This becomes equal to the Hubble time 
(3.08568 X 10^^sec)/(/i(l + z)^-^) at {1 + z) = 9/i° l Figure |l| shows the maximum 
allowed temperature^^ for an intergalactic medium heated impulsively as a function 
of the heating redshift. After being heated, the IGM cools by adiabatic expansion 
and Compton cooling. The 3 models shown in Figure |12| are an open model with 
Ho = 100 km s"^ Mpc"S fi^/i^ = 0.0125, and n = 0.2; a fiat model with = 50, 
= 1 and the same fis; and a "BARYONIC" model with Ho = 50 and [1 = ^3 = 0.2 
which has four times more baryons than BBNS models predict. For this increased 
value of the baryon abundance, the allowed kTmax is considerably less than the 1 keV 
suggested by the 300 km/sec velocity dispersion in galaxies. 

At redshifts z > Zy = 10^, there will be enough electron scattering to force the 
photons into a thermal distribution with a distortion instead of a ?/ distortion. 
However, the normal form for writing a /i distortion does not preserve the photon 
number density. Thus we should combine the /i distortion with a temperature change 
to give an effect that preserves photon number. The photon number density change 
with /X is given by 



1.5(1 + ni/ne)meC^ _ 7.4 x 10^^ f 



sec 



(21) 
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A similar calculation for the ener; 



■g; 



;y density shows that 



f/cx6(C(4)-/iC(3) + ...)- 



(23) 



In order to maintain = const, the temperature of the photon field changes by an 
amount AT/T = fj,({2)/{3C{3)). Therefore, the energy density change at constant N 
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Figure 12: The maximum allowed temperature for an IGM as a function of the heating 
redshift for two models with the BBNS baryon abundance (^2^/1^ = 0.0125) and one 
model with 4 times more baryons {Ho = 50, Q = Qb = 0.2). 
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AU /4C(2) C(3)^ 



U 



fi = 0.714/i. 



(24) 



V3C(3) C(4). 

Thus the FIRAS limit < 3.3 x 10"^ implies AU/U < 2.4 x 10"^. The "improved" 
form of the fi distortion, with a AT added to keep N constant, can be given as a 
frequency dependent brightness temperature: 



= To 1 + /i 



C(2) 
3C(3) 



+ 



(25) 



In this form it is clear that a /i distortion has a deficit of low energy photons and 
a surplus of high energy photons with respect to a blackbody. In this it is like the 
y distortion, but the crossover frequency is lower. The "improved" form of the fi 
distortion is plotted in Figure 



Finally, at high enough redshift the process of double photon Compton scattering 
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becomes fast enough to produce the extra photons needed to convert a distorted 
spectrum into a blackbody. Whenever a photon with frequency v scatters off an 
electron, there is an impulse oc hvjc transferred to the electron. This corresponds to 
an acceleration a oc hiP' /{niec) for a time interval At oc l/v. The energy radiated 
in new photons is thus oc e^h^v^ /{mic') which is oc ahv{hv / {meC^)Y ■ Since the 
rate of scatterings per photon is oc + zY ■, the overall rate of new photon 

creation is oc + z)^ while the Hubble time is oc (1 + z)^^. Thus the photon 

creation rate per Hubble time is oc (1 + 2;)^ and at earlier enough times a blackbody 
spectrum is produced no matter how much energy is transferred to the photon field. 
The emissivity for double photon Compton scattering has approximately the same 
spectral shape as free-free emission, since both arc due to impulsive accelerations. 
Thus the addition of photons to the radiation field can be described by 

9n 1 - e"^ / 1 \ 

^ = A—::^[-^-n] (26) 



dy Ve^ — 1 

where A is the ratio of photon creation via double photon Compton scattering or free- 
free emission to the increase of y. For double photon Compton scattering A scales like 
(1 + 2;). Because the photons are created primarily at low frequencies, n{x) approaches 
a blackbody for a; — 0, while the frequency shifts due to Compton scattering are 
simultaneously transferring photons to larger values of x. A is very small at redshifts 
where distortions can survive, so the photon creation occurs primarily at x <^ 1. For 
any n — (exp(a; -|- ^J^{x)) — 1)~^, 

so when /i(x) = /iexp(— Xo/x) and x,Xo,n and A are all <^ 1, implying n ~ n -|- 1 ~ 



This cancels the dn/dy due to photon addition, giving a quasi- equilibrium solution, 
if Xo — VA- With this form for n one finds a net photon addition rate of 

dN r odn , 

X -—ax 



dy J dy 



I 



2 , 1 - e- V 1 y ^ 

X A — n] ax 

^6 y^x _ 1 



= Afi/xo = n\fA. (29) 

Since the deficit of photons associated with is [I'k'^ j?)^ one finds a thermalization 
rate per unit y of 

■^In/X 3^1 nr—— 

— — = — — oc V 1 + (30) 
oy 
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Since {l+z)dy/dz oc (^5/1^(1+2;)^, the overall rate for eliminating a ^ distortion scales 
like ^Ish'^il + -z)^/^ per Hubble time. A proper consideration^*^ of this interaction of 
the photon creation process with the Kompaneets equation shows that the redshift 
from which 1/e of an initial distortion can survive is 



4.24 X 10 

which is zth = 2.3 x 10^ for the BBNS value of ^Bh^. 



^th = ,^ ,0^0.4 (31) 



8. Summary 

So far COBE has been a remarkably successful space experiment with dramatic 
observational consequences for cosmology, and the DIRBE determination of the 
cosmic infrared background is yet to come. The very tight limits on deviations of 
the spectrum from a blackbody rule out many non-gravitational models for structure 
formation, while the amplitude of the AT discovered by the COBE DMR implies 
a magnitude of gravitational forces in the Universe that is sufficient to produce the 
observed clustering of galaxies, but only if the Universe is dominated by dark matter. 
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